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We present a solution of the three-fermion problem in a harmonic potential across a Feshbach 
resonance. We compare the spectrum with that of the two-body problem and show that it is 
energetically unfavorable for the three fermions to occupy one lattice site rather than two. We also 
demonstrate the existence of an energy level crossing in the ground state with a symmetry change of 
its wave function, suggesting the possibility of a phase transition for the corresponding many-body 
case. 



I. INTRODUCTION 



II. METHODS 



Ultracold atoms, tuned with Feshbach resonance, offer 
a great opportunity to study strongly correlated many- 
body physics in a controlled fashion. For such strongly 
interacting systems, in general there is no well-controlled 
approximation method to solve the many-body physics. 
Exact solution of few-body problems plays an important 
role in understanding the corresponding many-body sys- 
tems. Few-body (three- or four-body) problems have 
been solved for strongly interacting bosons or fermions 
in free space jl|, in a quasi-one-dimensional configura- 
tion (2j, and in a three-dimensional (3D) harmonic trap 
for bosons Q or fermions in the unitary limit 0, Q . 

In this paper, we add to this sort of examples by ex- 
actly solving the three-body problem for strongly inter- 
acting two-component fermions in a 3D harmonic trap 
across resonance. This work has two main motivations: 
Firstly, the situation considered here is relevant for ex- 
periments where one loads strongly interacting two com- 
ponent fermions into a deep 3D optical lattice @. For 
each site that can be approximated with a harmonic po- 
tential, one could have two identical fermions (spin-|) 
strongly interacting with another distinct fermion (spin- 
|). The three-body problem for equal mass fermions 
turns out to be very different from the corresponding case 
for bosons. Instead of a hierarchy of bound Efimov states 
for bosons we show that there is always a significant 
energy penalty for three fermions to occupy the same 
lattice site (Ttl) instead of two (|J. + |). This result 
justifies an important assumption made in the derivation 
of an effective many-body Hamiltonian for this system 
[3]. Secondly, we analyze the ground-state structure of 
the three-body problem and show that as one scans the 
3D scattering length, there is a level crossing between 
the lowest-lying three-fermion energy eigenstates, which 
have s- or p-wavc symmetries respectively in the limit 
as two atoms are contracted to form a dimer. This level 
crossing with a symmetry change may correspond a quan- 
tum phase transition in the many-body case where one 
has spin-polarized fermi gas loaded into an optical lat- 
tice. The latter system with polarized fermions has raised 
strong interest recently in both theory and experiments 

aa. 



The method here is based on manipulation of a 
Lippmann-Schwinger equation for the wavefunction, with 
the formalism similar to the one presented in Ref. (2j. 
As is standard, we separate the center-of-mass degree of 
freedom via an orthogonal transformation of variables, 
leaving the trapping potential diagonal in the new coordi- 
nates. The Schrodingcr equation for the relative degrees 
of freedom is then 
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where too is the atomic mass, lo is the trap frequency, 
y is the vector between the two ] fermions, \/3x/2 is a 
vector from the center of mass of the two ] fermions to 
the J. fermion, and r± = \/3x/2 ± y/2 are the vectors 
from the | fermion to each of the two | fermions. We ap- 
proximate the short-range interaction between fermions 
with the usual zero-range pseudopotential V (r) = 
4 ^ a S (r) (r-), where a is the 3D s-wave scattering 
length tunable through the Feshbach resonance. 

The above contact interaction is equivalent to impos- 
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ing boundary conditions 'F (x, y) ~ =p 4 

for r± — > 0, where the rj_.± = x/2 =F \/3y/2 are pro- 
portional to the distances between the center of mass of 
an f J, pair and an f fcrmion. The overall =F sign ensures 
the antisymmetry of the wavefunction upon swapping the 
identical fermions. The undetermined function / (tj_,±), 
after a rescaling of the argument, is the relative atom- 
dimer wavefunction that results when two of the fermions 
form a tightly bound pair. Solving for this asymptotic 
wavefunction / (r± y ±) fully determines 'F (x, y). 

Since V (r±) only acts at r± = 0, we use the asymp- 
totic form of ^ when computing V (r±) 'F. The formal 
solution can then be written as 
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where the two-particle Green's function is given by 



Gf (x,y;x',y')=£ 



^ (x)^ A2 (y)^ (x')^ a (y') 
£ Al +E X . 2 -E 

(3) 

The ^/i Ai (x) are the single-particle eigenfunctions with 
eigenenergies for the reduced mass mo/2. Here we 
use spherical coordinates for the three-dimensional har- 
monic trap, so the quantum numbers are A = (n,l,m), 
n = 0,1,2,...; I = 0,1,2,...; m = -I, -I + 1, I - 1,1. 
The eigenenergies are E\ = (2n + I + 3/2) Too and the 
eigenfunctions are ip\ (r) = R n i (r) Y™ (9, 4>), where the 
Y, m (9, <fi) are the standard spherical harmonics. The ra- 



dial wavefunction is given by [11[ 
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is the length scale of the trap and the 



L\ (r) are associated Laguerre polynomials. 

We can make use of the invariance of G E and the 
integration measure under an orthogonal transformation 
of variables to rewrite Eq. in terms of r = r_ and 
r 1 = r 1 ._, 



*(r,r ± ) = — / dr' ± f(r' ± 
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We can also decompose the asymptotic atom-dimer wave- 
function in terms of the complete set of single-particle 
wavefunctions, f (r±) = J2x fx^Px ( r J-)- Then Eq. ([5]) 
becomes 
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FIG. 1: Energy vs. inverse scattering length. 
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FIG. 2: Difference between three-fermion energy and two- 
fermion energy plus one-fermion energy vs. inverse scattering 
length. 



work, we obtain 
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where Ge is the single-particle Green's function 
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fl2j |. and t/ is the confluent hypergeometric function. 
Note that the three-fermion wavefunction is fully deter- 
mined by / (rj_) and that if we consider Eq. ([6]) in the 
limit as r — > 0, we obtain a self-consistent equation for 
/( r J_) by using the boundary conditions. After some 
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III. RESULTS 



We anticipate that the low-energy physics should be 
contained in a truncated Hilbcrt space containing only 
the lowest few asymptotic atom-dimer energy levels. 
Then Eq. ([8]) is easily solved numerically. We have 
checked that, indeed, the solution for the ground state 
and the first excited manifold become insensitive to the 
cutoff, as long as the first four or five atom-dimer energy 
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FIG. 3: (Color online) Relative coordinates for the three- 
fermion problem. In terms of variables used in Eq. (O, r — 
|r|, R= V3\r±\/2. 




(a)n = 0, I = 0, m = 




(b)ra = 0, I = 1, m = ±1 



(c)n = 0, I = 1, m. = 

FIG. 4: Contour plots of r 2 |* (r, i£, 0)|*_,. o as a function of R 
and for r/d ~ and d/a ~ 0. 



levels are included. For all results presented in this pa- 
per, we have kept the first five energy levels. Because of 
the degeneracy of the excited levels, Eq. ([8]) becomes a 
35 x 35 matrix equation. For a given energy, we solve nu- 
merically to obtain the corresponding scattering length 
and eigenstate. By sweeping through a range of energies, 
we map out the spectrum shown in Fig. [TJ At unitarity, 
our result for the low-energy spectrum agrees with the 
analytic result of Ref. Upon careful inspection, one 
can discern the presence of level crossings. More usefully, 
in Fig. [2] we display the difference between the energy of 
three fermions in a single lattice site and the energy of 




(a)n = 0, I = 0, m = 




(b)n = 0, I = 1, m = ±1 




(c)n = 0, I = 1, rn = 



FIG. 5: Contour plots of r 2 R 2 |sin6>* (r, R, 8)\ 2 as a function 
of R and 6 for r/d = 2 and d/a ~ 0. 

two fermions (||) in the site and the extra fermion alone 
in a separate site. Here we have used the well-known 
exact solution for the two-body energy [HJ , 
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Figure [5] is our main result. There are two main fea- 
tures we would like to point out. First, it is clear that it 
is energetically favorable for atoms in a lattice to arrange 
themselves such that there are less than three atoms per 
site, regardless of the scattering length. This has already 
been assumed in the derivation of an effective many-body 
Hamiltonian for atoms in an optical lattice across a Fes- 
hbach resonance @, and is confirmed by Fig. Second, 
the level crossing in the ground state is now quite evi- 
dent. On the positive scattering length side of the cross- 
ing (the BEC side, where the many-body system forms a 
Bose-Einstein condensate of bound dimers) , the ground 
state is nondegenerate. On the other side (the BCS side, 
where the many-body system forms a Bardeen-Cooper- 
Schricffcr superfluid of atomic Cooper pairs), it is triply 
degenerate. Other crossings appear in the excited spec- 
trum, although to obtain quantitatively accurate results 



4 






FIG. 6: (Color online) Normalized probability density distribution functions of variables defined in Fig. [3]for various scattering 
lengths. The solid line is for the n = 0, I = 0, m = state; the dashed line is for the n = 0, I = 1, m = ±1 state; and the 
dotted line is for the n = 0, 1 = 1, m = state. 



for these one should include higher modes when solving 
Eq. ©. 

The origin of the level crossing is the differing symme- 
tries of the eigenstates. For very small, positive scatter- 
ing length (deep BEC side), formation of tightly bound 
dimers is favorable, so the state should behave as the 
ground state of the relative atom-dimer motion, which 
has s-wave symmetry. For very small, negative scatter- 
ing length (deep BCS side), the atoms are essentially 
non-interacting, so the ground state comprises two atoms 
(ti) in the ground state of the trap plus the third in the 
first excited state of the trap (which is triply degener- 
ate) because of Pauli exclusion. So, on the deep BCS 
side, the ground state has p-wave symmetry. Due to the 
rotational symmetry of a spherical harmonic trap, the 
total angular momentum of the three particles should be 
a conserved quantity. However, from the above analysis, 
this quantity has different values for the ground state in 



the deep BEC and deep BCS limits. Therefore, there 
must be a ground-state level crossing for this system as 
one scans the scattering length. If one considers multiple 
lattice sites with each site having on average two spin 
t and one spin j atoms (which could be realized with 
polarized fermions in an optical lattice with appropriate 
filling number and population imbalance), as the three- 
body problem has a level crossing with different ground 
state degeneracies in the BCS and the BEC limits, there 
could be a corresponding quantum phase transition for 
this many-body system (with small tunneling between 
lattice sites) as one scans the scattering length. 

The wavefunction given in Eq, @ does not generally 
have definite relative angular momentum for any two 
fermions. However, in the limit as the distance between 
two distinguishable fermions goes to zero, the wavefunc- 
tion takes on the symmetry of the asymptotic atom-dimer 
wavefunction in the remaining coordinates. This is a rel- 
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ative angular momentum cigenstate due to the spherical 
symmetry of the limiting case. With the relative coordi- 
nates defined by Fig. [3l the symmetry of the wavefunc- 
tion as a function of R and 9 in the limit as r goes to zero 
is shown by Fig. [4] For finite r, the m = wavefunc- 
tions are affected by the asymmetry and take nontrivial 
shapes. We have plotted an example in Fig. Note 
that in this figure we include a factor of sin 9 since the 
wavefunction itself diverges at R = r /2, 9 = tt according 
to the boundary conditions we have imposed. On the 
BEC side, the lobes are tightly bunched near 9 = tt, but 
on the BCS side they spread around as expected. 

In general, the eigenstates are rather difficult to vi- 
sualize, since they depend nontrivially on three spatial 
variables as well as the scattering length. However, one 
can get some idea of the evolution of the eigenstates from 
Fig. [6l which shows the normalized probability density 
as a function of the variables introduced in Fig. [3]for var- 
ious scattering lengths. In each subplot we have numeri- 
cally integrated over the other two variables to obtain a 
one-dimensional function. The dimer size, r, clearly de- 
creases in size as one enters the BEC regime, developing 
a strong peak at the origin. However, as is clear from 
Fig. [3l there are two ways to form a tightly bound dimer 
(due to the two identical spin f atoms) and we have arbi- 
trarily chosen one to define the origin r = 0. So it is not 
surprising that we see a more diffuse second peak at large 



distance r, corresponding to the dimer forming between 
the spin J, and the other spin | atom. This is also the 
meaning of the spike at 9 — tt on the BEC side. 



IV. SUMMARY 

We have found the low-lying energy levels of three 
fcrmions in a harmonic trap and examined the corre- 
sponding wavefunctions. The ground state has s-wave 
symmetry on the BEC side of Feshbach resonance and 
has p-wave symmetry on the BCS side. In the resonance 
region there is a level crossing, which may indicate a 
phase transition in the corresponding many-body case. 
We also note that, in the vicinity of resonance, the en- 
ergy of three atoms in a single site is greater than their 
energy if they are in two sites, with a gap on the order 
of the trap spacing, validating the approximation in Ref 

0- 
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Note added: After com plet ion of this work, we became 
aware of a recent work |13[ which treats the trapped 
three-fcrmion problem with a different approximation 
method. 
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